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Abstract 

We study renormalization group running effects on neutrino mixing patterns when a 
(type I) seesaw model is implemented by suitable flavour symmetries. We are partic- 
ularly interested in mass-independent mixing patterns to which the widely studied 
tribimaximal mixing pattern belongs. In this class of flavour models, the running 
contribution from neutrino Yukawa coupling, which is generally dominant at ener- 
gies above the seesaw threshold, can be absorbed by a small shift on neutrino mass 
eigenvalues leaving mixing angles unchanged. Consequently, in the whole running 
energy range, the change in mixing angles is due to the contribution coming from 
charged lepton sector. Subsequently, we analyze in detail these effects in an explicit 
flavour model for tribimaximal neutrino mixing based on an discrete symmetry 
group. We find that for normally ordered light neutrinos, the tribimaximal prediction 
is essentially stable under renormalization group evolution. On the other hand, in 
the case of inverted hierarchy, the deviation of the solar angle from its TB value can 
be large depending on mass degeneracy. 
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1 Introduction 



The seesaw mechanism [1] explains the lightness of left-handed (LH) neutrinos in a simple 
and elegant way. However, the seesaw models usually contain many more parameters 
than low energy observables. In order to economically describe the neutrino oscillations, 
new ingredients are needed: the introduction of a horizontal symmetry could improve the 
situation and indeed many examples of this kind are present in literature. A common 
feature of these models is to provide a description of neutrino masses and mixings only 
at a very high energy scale. On the other hand, for a comparison of experimental results 
with the high energy predictions from flavour symmetries, it is important to evolve the 
observables to low energies through renormalization group (RG) running. In general the 
deviations from high energy values due to the RG running consist in minor corrections, 
but the future improvements of neutrino experiments could hopefully bring the precision 
down to these small quantities. 

In the MSSM context, the RG effect on the light neutrino mass operator m u in the 
leading Log approximation can approximately be parametrized as 

TTlv (lower Energy) IjjJe (higher Energy) Ju Je 

where Ijj is a universal contribution, J e , that is proportional to Y^Y e with Y e the charged 
Yukawa coupling, is always flavour dependent and, J v may depend on flavour or not. 
The running contributions to m v from J e are similar to those below the lightest right- 
handed (RH) neutrino mass and it is as known well under control. In fact large RG 
effects can be expected only in the case of a degenerate light neutrino spectrum or for 
large values of tan/? [2]. The running contribution from J u , however, is generally more 
complicated depending non-trivially on the neutrino Yukawa coupling Y v . Moreover, the 
last contribution can even be the dominant one since Y v is expected to be of order one [3] . 
Similar conclusions can be driven in the SM context, when assuming a flavour origin of 
the neutrino mass term m„. In seesaw models implemented by flavour symmetries, Y v is 
usually subject to constraints in order to efficiently describe the observed neutrino mixing 
structure. We should expect that these constraints have also some important impacts 
on running effects. In this paper, we will show that, under quite general assumptions, 
flavour symmetries imply a J v contribution which has not effects on mixing angles. As a 
consequence even including energy ranges above and between the RH neutrino scales, the 
contribution from the charged lepton Yukawa coupling J e on mixing patterns dominates. 

In the first part of the paper, we will describe, in a very general context, two kinds of 
interesting constraints on Y u from flavour symmetries and then analyze their impact on 
running effects. We first consider flavour models in which Y v is proportional to a unitary 
matrix. It is the case, for example, when the RH singlet neutrinos or the charged leptons 
are in a irreducible representation of the flavour group Gf. 

Then we extend this constraint to a more general class of flavour models in which the 
mixing textures are independent from the mass eigenstates. As a general result, we find 
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that in this class of models, the effect of J v can always be absorbed by a small shift on 
neutrino mass eigenvalues and the mixing angles remain unchanged. This conclusion is, 
in particular, independent both from the specific mixing pattern implied by the flavour 
symmetry and the basis where we are working. 

Mass-independent mixing textures usually exhibit a underlying discrete symmetry na- 
ture. For example, the tribimaximal (TB) mixing pattern [4] which is defined by 



belongs to this class of mixing textures. TB mixing pattern, widely studied in the last 
years, provides a very simple first order description of the existing oscillation data [5-7]. 
In fact, it has been realized that the TB mixing matrix of eq. (1) can naturally arise as 
the result of a particular vacuum alignment of scalars that break spontaneously certain 
discrete flavour symmetries. A class of economical flavour models which can naturally 
explains TB mixing is based on A4 discrete symmetry [8] [9] [10,11]. Subsequently the 
symmetry group A 4 has also been extended to the group T' [12] and S4 [13, 14] to cover 
a reasonable description also for quarks. Another important historical example of mass- 
independent mixing scheme is the bimaximal (BM) mixing (see [15] for a recent revival of 
BM in the context of the discrete symmetry S4 and for an up-to-date list of references). 
The running effects on BM [16] and TB [17] mixing patterns have already been studied in 
literature without, however, consider an explicit realization based on flavour symmetry. 

Then, as an explicit example, in the second part of the paper, we describe in detail the 
RG effects on the TB mixing texture in the lepton flavour model proposed by Altarelli- 
Feruglio (AF) [8]. At leading order, this model contains less parameters than the case of 
a general TB pattern considered in [17]. Then our result may not a priori reproduce the 
same results obtain in [17]. In our analysis, the RG corrections to mixing angles and the 
phases are discussed as functions of the lightest neutrino mass and the type of spectrum. 
The analysis has been performed both in the Standard Model (SM) framework and the its 
minimal SUSY extension (MSSM). 

The outline of the paper is organized as follows. In section 2, we analytically discuss 
the RG effects on the neutrino mass operator m v in a type I seesaw model. In section 
3, we characterize two general classes of flavour models in which J v has no effects on 
mixing angles. Then we turn to RG effects from J e on TB mixing pattern. In section 4 
we introduce the AF model and show its main features including also the next-to-leading 
(NLO) contributions coming from higher dimensional operators. In section 5, a more 
detailed numerical analysis on the impact of RG running in the AF model is given. The 
result has been compared also with the NLO contributions. In the end in section 6 we 
conclude summarizing our main results. 




^fljl l/y/3 

-l/Ve i/Vs -1/V2 

-l/VE l/y/3 +1/V2 



(1) 
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2 RG effects on neutrino mass operator m u 

In this section we begin to analyze, in a general context, the RG equations for neutrino 
masses below and above the seesaw threshold, both in the SM and in MSSM extended with 
three right-handed neutrinos. For definiteness, we consider the following lagrangian 



e cT Y e HH + v cT Y v EH + v 



cT 



M R v c + h.c. 



(2) 



where I are the LH lepton doublets, e c the RH charged lepton singlets and H(H = ia 2 H*) 
is the Higgs doublet. In the supersymmetric case, this lagrangian should be identified to a 
superpotential where H(H) is replaced by h d (h u ) and all the fields are instead supermul- 
tiplets. In what follows we concentrate only on the SM particles, postponing the study of 
their supersymmetric partners elsewhere: for this reason in our notation a chiral superfield 
and its -R-parity even component are denoted by the same letter. 

Given the heavy Majorana and the Dirac neutrino mass matrices, M R and rri£> = Y v v/y2 
respectively, the light neutrino one is obtained from block-diagonalising the complete 6x6 
neutrino mass matrix: 

„2 _ 

(3) 



m. 



—Y T M~ X Y 
2 1 v m R ' 



where v refers to the VEV of the Higgs field, (H) = v/^/2 (v ~ 246 GeV). The equivalent 
relation in the supersymmetric case is achieved by replacing v with v u . In our notation 
the VEV of the Higgs fields h u and h d are given by {h Ujd ) = v u>d /y/2, with \Jv\-\-v\ = v. 
The matrix m v is modified by quantum corrections according to the renormalization group 
equations (RGEs) widely studied in literature [3]. For completeness, in Appendix A, we 
report the full RGEs for all the interested quantities in the running. Here, we explicitly 
give RGEs only for the Yukawa couplings Y e , Y v and the RH neutrino mass Mr, relevant 
for the analysis on this section: in the SM context they are written as 



3YjF u + W}Y d + YjY u + Y}Y e 



16tt 2 4y = Ye{-Y}Y e - -Y}Y U + Tr 
dt [2 e 2 u 

167r 2 -Y, = Y v \ -YjY u - -Y}Y e + Tr \3YjY u + 3V]V d + Y$Y V + Y}Y e 
16vr 2 ^M fl = (Y U YJ) M R + M R {Y u Y}f , 

while in the MSSM they are slightly different, 



4? x 4^ 2 



20- 



4^ 



(4) 



Y e I 3YjY e + YjY u + Tr 3Y H % + YjY t 



16tt 2 4y 
dt 

l6n^Y u 
dt 

16n 2 j t M R = 2(Y V YJ) M R + 2M R (Y„Y}f 



Y v 3YJY U + Y}Y e + Tr 3YjY w + Y$Y V 
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:gl - *gl 



-gl - 3gl 



(5) 
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where t := ln(/x//i ) and Y u ^) is the Yukawa coupling for the up- (down-) quarks. 1 From 
these RGEs, one can obtain the RG evolution for the running composite operator m v {[j) 
defined in eq. (3). In order to analytically study the change of m v {ix) from high to low 
energy, it is useful to work in the basis in which the Majorana neutrino mass is diagonal 
and real, Mr = diag(Ms, Mm, Ml). The mass eigenvalues can be ordered as Ms < Mm < 
M^. Furthermore, we can divide the RG effects in three distinct energy ranges: from 
the cutoff A of the theory down to Ml, the mass of the heaviest RH neutrino; from Ml 
down to Ms, the mass of the lightest RH neutrino; below Ms down to A, which can be 
either rriz, considered as the electroweak scale, or msusY, the average energy scale for the 
supersymmetric particles. 

Af — ► Ml- Above the highest seesaw scale the three RH neutrinos are all active and the 
dependence of the effective light neutrino mass matrix from the renormalization scale 
\x is given by mean of the \i— dependence of Y v and M R : 



mM = --Yj^M^YM. 



(6) 



Then from the RGEs in eqs. (4, 5), it is not difficult to see that the evolution of the 
composite operator m v is given by: 



16?^ 



m u + m u (c e Y}Y e + C U Y^ + 



a m. 



(7) 



with 



and 



«5M 



&MSSM 



2 

C e = C u = 1 
= 2Tr 



a 



in the SM 
in the MSSM 



2Tr 



3YJY U + 3YjY d + YJY V + Y^Y e 



3Y^Y U + YJY U -\g\-Qgl. 
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10 9l 2 92 



(9) 



Ml — ► Ms- The effective neutrino mass matrix m v below the highest seesaw scale can 
be obtained by sequentially integrating out z/£ with n = L,M, S: 



v ( in) < n V (n) i (n) 
m v = -— rn + 2YjM- l Y v 



(10) 



where k is the coefficient of the effective neutrino mass operator HHHH. From the 
(tree-level) matching condition, it is given by 



'IF 



the GUT normalization, such that g 2 = g and gi = Uhj2>g' . 
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(L) 

which is imposed at fi = M n . At Ml, the 2x3 Yukawa matrix Y u is obtained by 

simply removing the L-th row of Y v and the 2x2 mass matrix Mr is found from 
Mr by removing the L-th row and L-th column. Further decreasing the energy scale 

(M) (L) 

down to Mm, Y v is a single-row matrix, obtained by removing the M-th row from Y u , 

(M) 

and Mr consists of a single parameter, found by removing the M-th row and M-th 

(L) (S) (S) 

column from Mr. Finally at Ms, Y v and Mr are vanishing. 

In the SM, the two parts which define m v in eq. (10) evolve in different ways. We 
can summarize the corresponding RGEs as follows: 



(n) 

dX 



dt 



1 (n) (n) 



16vr 2 — = ( -Y$Y V - -y e +y e ) X + X ( -Y},Y V - -Y]Y t 



(n) (n) / ]_(n) (n) 



e 1 e 



(n) (») 

otxX 



(12) 



where 



(n) 



(n) 



2Tr 
2Tr 



(«),(«) 



3y u ty + 3yjy, + yjy + y e ty e 



3y u ty + 3y;y, + yjy + y e +y 



3(7| + Ah 



10 ^ 2^ 2 



(13) 



with A# the Higgs self-coupling. 



(n) (n) (n) 



In MSSM the running of Ac and of Y^M^Y V is the same and therefore we can write 
2 dm v 



16tt" 



dt 



T 



(n) (n) \ 1 f (n) (») \ 

y e ty e + yjy m, + m v (y}y £ + yjy J + 



(n) 



where 



a - 



2Tr 



(«),(») 



6 2 R 2 

g01 - • 



(14) 



(15) 



Ms — > A. For energy range below the mass scale of the lightest RH neutrino, all the i/° 

(S) (S) 

are integrated out and Y v and Mr vanish. In the right-hand side of eq. (10) only the 
term ac is not vanishing and in this case the composite operator m u evolves as: 



with 



167T 2 ^ = (c e Y}Y e ) T m u + m v (c e Y}Y e ) + 



(S) 

a m v 



(S) 

®SM 



(S) 

OtMSSM 



2Tr 
6Tr 



3yjy„ + 3y;y d + y e ty e ] - 3<7 2 + a^ 

6 



YlY, 



- t91 - ®gi 



(16) 



(17) 



2 We use the convention that the Higgs self-interaction term in the Lagrangian is —\h{H^H) 2 /A. 
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2.1 Analytical approximation to RG evolution of 

Now we analytically solve the RG equations for m v in the leading Log approximation. All 
the Yukawa couplings Y}Yi for % — u, e, u, d are valuated at their initial value at the cutoff 
Af. Furthermore we will keep only the leading contributions from each Y}Yi term, for 
i = e,u,d, i.e. \y T \ 2 , \yt\ 2 and \yb\ 2 respectively. The RG corrections to these quantities 
would contribute to the final result as sub-leading effects and we can safely neglect them 
in the analytical estimate. 

In the MSSM context, the general solution to eqs. (7), (14) and (16) have all the same 
structure, which is approximately given by 



ni 



v (lower Energy) 



~ 1 U J e v '"^(higher Encrg; 



;is) 



where Iy, J e and J u are all exponentials of integrals containing loop suppressing factors 
and as a result they are close to 1. Note that ly is a universal contribution defined as 



exp 



16tt 



a at 



(19) 



where the integral runs between two subsequent energy scales and we have extended the 

(n) (A) 

definition of a by identifying a = a in order to include the range from A down to Ml- J e is 
the contribution from charged lepton Yukawa couplings which is always flavour- dependent 
and is given by 3 

1 



Je 



cxp 



167T 2 



Y]Y P dt 



(20) 



Finally, J v is the contribution from the neutrino Yukawa coupling 



J u = exp 



1 f (n).(n) 

— ; / yJy u dt 

16tt 2 / v 



(21) 



(n) (A) 

where also here we have extended the definition of Y v by identifying Y v with Y v in order to 
include the range between A and Ml. Differently from J e , J u can be flavour-dependent or 
not. 

In the SM context, the RG effect does not factorize, due to the different RG evolution 

of k and YJ M R Y u between the seesaw mass thresholds. However eq. (18) applies also 
to the SM context when m v is a result of a flavour symmetry: in this case, by a suitable 

redefinition of the mass eigenvalues, the sum k + YjM^ Y y after the RG evolution has 
exactly the same flavour structure of m u (higher Energy) • For the purposes of the present 
discussion we simply assume that eq. (18) is valid also in the SM context and an explicit 
example will be proposed in section 3.3. 



3 In eq. (20), the combination Y^Y e should enter with Y e instead of Y e , as one can see from the RGEs 
in Appendix A. In our approximation, however, they coincide. 
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Expanded J e and J u in Taylor series and summing up (18) on several energy ranges 
one can approximately calculate the neutrino mass at low energy as 



m v{ x) ~ I v (m v{A ) + Am< Je) + Am^ 1 ) 



(22) 



where the low energy scale A is rriz in the case of SM and mgusY for MSSM. The explicit 
form of the universal part Iy is given by: 



rSM 



1 x exp 



16tt 2 



(-^-^ + %,| 2 )ln^ + (^ + ^ + A„).„ 



Mc 



m z 



/ 2 ln^ + 41n^ + 71n A ' 



M 



M 



Mr 



rMSSM 
l U 



1 x exp 



16tt 2 



:gj-6gl + 6\y t \ 2 1 hi 



A 



(23) 



/ 



m SUSY 



+ 



+ y 2 2 In 



M 



Mc 



+ 4 In 



Mr 



M 



+ 8 In 



A/ 

Mr 



Ami J ^ is the the contribution from J e and can easily be calculated as: 



Am™ = m KA ) diag(0, 0, A r ) + diag(0, 0, A T )m^ 



(A) 



(24) 



where the small parameter A T is given by 
A r = 



3m 2 : , A 

tt-t: m ■ 



lQn 2 v 2 mz 



m A 
A r = _^_(l + tan 2 /3)ln- 



8tt 2 v 



in the SM 



in the MSSM 



(25) 



fflSUSY 



where tan (3 is the ratio between the VEVs of the neutral spin zero components of h u and 
hd, the two doublets responsible for electroweak symmetry breaking in the MSSM. On the 
other hand, the contribution from J u , Ami Jv \ non trivially depends on the neutrino Yukawa 
coupling Y v which cannot be determined by low energy observables without additional 
ingredients. In section 3, we will analyze strong impacts of the flavour symmetries on J v . 
But before proceeding, we can first give a naive estimate of various running contributions 
to neutrino mass. 



2.2 Naive numerical estimate 

Here we numerically estimate the contributions of running, encoded in Am„, above the 
seesaw scales (Am !y ( high )), those between the seesaw scales (Am I/ ( seesaw )) and below the 
seesaw scales (Am v n ow )) on the light neutrino mass matrix. One should expect that the 
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mixing angles will be affected by the same quantities. We assume that flavour symmetries 
have not effects on Y v and then, without cancellations, 

(«),(«) 

YX ~ YjY v = Oil) . 

As a result, both in the energy range above and between the seesaw scales, the running 
contributions from Y u which are all encoded in J u would dominate over the one from Y e , 
encoded in J e . Summing up all the two contributions we have 

|Am v (high)| + |Am i , (sccsaw) | ~ |Am^| ~ — ^ In (26) 

107T IVlg 

both in the SM and in the MSSM contexts. Below the mass of the lightest RH neutrino 
Mg the contribution J v is absent and we should roughly obtain 

3m 2 M s . 

Arrwiow) ~ , n o o m m the SM 

167i 2 v 2 m z (97 , 

|Attw 1ow) | ~ — ^(l + tan 2 /3)ln S — in the MSSM . 

Now we can compare the two contributions Ami J ^ and Am I/ (i ow ). In general one assumes 
a seesaw scale between 10 12 GeV and 10 14 GeV and a cutoff scale approximate to the 
grand unified scale Af ~ 10 16 GeV. For the numerical estimate we use Ms ~ 10 12 GeV, 
msvsY ~ 10 3 GeV and mz ~ 100 GeV. Then we obtain 4 

~ 2.5 x 10" 4 in the SM 

(28) 

[1 + tan 2 /5) ~ 1.5 x 10" 4 tan 2 /3 in the MSSM . 

We can conclude that both in the SM and the MSSM, the contribution from J v always 
dominates over Am^^) even for large tan/5 (we consider tan/3 = 60 as the maximal 
value) . 

The previous naive estimate for RG effects on the neutrino mass operator, however, 
cannot always be transfered to the change in mixing angles A^- for a particular model 
building. Indeed, as we will see in a moment, quite frequently flavour symmetries imply 
a J v which is flavour-independent or has no effects on mixing angles. If it is the case, 
the situation changes drastically: even if Am y ( low ) can still be dominated by Am I/ ( seesaw ) , 
A0y(iow) should dominates over A%( seesaw ). 



Am^iow) 


3m 2 8 

rsj 


|Amt Jl/) | 


v 2 5 


Am^iow) 


2m 2 7 

rsj 


\Aml Jv) \ 


v 2 5 



3 Flavour symmetries and RGE effects 

In the present section, we will apply the general results of the RG evolution of the neutrino 
mass operator m u to models beyond the Standard Model, where a flavour symmetry is 

4 We take the pole mass for the lepton r [18]: m T = (1776.84 ± 0.17) MeV. 
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added to the gauge group of the SM. The main aim is to track some interesting connections 
between the running effects and how the flavour symmetry is realized in nature. 

In a given basis, Y^Y e and m v can be diagonalised by unitary matrices, U e and U u , 
respectively The lepton mixing matrix is given by C/pmns — U\U V . The analysis of 
how £/pMNS changes with the RG running has already extensively performed [2] in the 
context of SM and MSSM. On the other hand, only few studies [3] are present in literature 
considering the presence of additional RH neutrinos, which originate the type I seesaw 
mechanism. Here we develop a general RG analysis for seesaw models in which the lepton 
mixing matrix £7pmns is dictated by a flavour symmetry Gf. It is a common feature in 
flavour model building that Gf must be spontaneously broken in order to naturally describe 
fermion masses and mixings, as we will see in section 4. Here, we simply assume that Gf is 
spontaneously broken by a set of flavon fields $ at a very high scale. The symmetry group 
can be discrete or continuous, global or local (or even a combination of them). Suppose 
that, at leading order, the neutrino mixing matrix is given by Uq which differs from £/pmns 
by subleading contributions ~ ($)/Aj where Af is the cutoff scale of the flavour symmetry 
Gf. We will begin with some general assumptions on Uq without however specifying its 
form. Then we will move to specialize in a concrete case in which Uq is given by the TB 
mixing pattern. 

3.1 Running effects on neutrino mixing patterns 

As described in Sec. (2) the relevant running effects on m„ are encoded in the combinations 
YJY e and YjY v . Furthermore, we observe that a relevant contribution to the running of 
YjY e is encoded by YjY p . We perform the analysis in the basis in which the charged leptons 
are diagonal, then at high energy we have 

Y^Y e = dmg(mlm 2 m 2 T )^ . (29) 

From now on, we will use v in the notation of the SM and in order to convert similar 
expressions to the MSSM, it is sufficient to substitute v with v Ut d, when dealing with 
neutrinos or charged leptons, respectively. Naturally, this simple form should change when 
evolving down to low energies. The running effect of YJY e on m u is of second order and we 
can safely forget it. However it can generate a non trivial U e and consequently introduces 
additional corrections to £/pmns- We will return to this effect in section 3.2. 

Since flavour symmetries impose constraints on Y v , they should have some impacts 
also on running effects. In this paper we are interested in two classes of constraints. The 
first class is characterized by Y v proportional to a unitary matrix and in the second one we 
assume that m v can be exactly diagonalized by Uq according to 

rh u = UqUIvIIq (30) 

where rh u = diag(mi, 7712, ms) with rrii positive and Uq is a mass-independent mixing 
pattern enforced by the flavour symmetry Gf. Independently from the way Gf is broken, 
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we will show that, in this second case, the neutrino Yukawa coupling in the basis of diagonal 
RH Majorana neutrinos, which we indicate as Y v , has the following simple form 

Y v = iDUl (31) 

where D = diag(±V2m 1 M 1 , ±y / 2m 2 M 2 , ±y/2m 3 M 3 ) jv. 

3.1.1 YtY u = Y V Y\ ~ 1 

As already pointed out, YjY u plays an important role in determining the running effect 
on m v . In the case of a Y v proportional to a unitary matrix the study of RG evolutions 
becomes quite simple. YjY u ~ 1 or Y V Y^ ~ 1 is quite frequent in the presence of a 
flavour symmetry. It is, for example, a consequence of the first Schur's lemma when I or 
v c transforms in a irreducible representation of the group Gf [19]. 

Between the energy scales A/ and Ml one has J v oc 1 then its contribution to m v is 
global. Now we demonstrate the contribution from J v is universal also between Ml and 
Ms if YjY u ~ 1. We first consider the running effect on defined in eq. (11) which is 
proportional to 

(n) An) 

(YjY u ) k + symmetrization . 
In components, for the first term we have 

rriy^l k 

where we have used the unitary condition Ylk^X^mkiYv)^ = S mn . As we expected J v has 

(n) 

no non-global effects on k. Now we move to analyze the second term in (10). The running 
effect on this term is proportional to 

(n) (n) (n) (n) (n) 

(Y^Y V ) T Yj M^ l Y v + symmetrization, 
Similarly as before, using the unitary condition, we obtain in components 

E E(^ T )™(^)^( 1 ^ r )^( M ^" 1 )^(^)« = E W)*(^ _1 )«c^)«j . 

(n) (n) (n) 

which coincides with the starting observable, Yj M^ Y v . In the previous expression, the 
meaning of the apex (n) is understood. A similar result is valid for the symmetrization 
part. Then we conclude that m„ does not change the flavour structure under J u if Y v is 
proportional to a unitary matrix. This result in completely general and holds in any basis, 
not only in the hatted one. In this case, the only flavour-dependent RG contribution to 
m u is encoded in J e . 
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3.1.2 The general case Y„ = iDU^ 

Now, we consider a class of constraints on Y u weaker than the previous one but still very 
interesting for models based on flavour symmetries. We simply assume that at high energy 
m v can exactly be diagonalized by Uq and Uq is a mass-independent mixing pattern. The 
TB mixing pattern, independently from what is the underling flavour symmetry, is one 
of the examples of this class of models. Other examples are given by flavour symmetries 
which give rise, at leading order, to the Bi-maximal mixing pattern [15], to the golden ratio 
mixing [20] and some (but not all) cases of the trimaximal mixing [21]. 

In the base where the RH Majorana mass Mr is diagonal eq. (30) can be written as 

ro„ = -y UfYj M~ l %Uq ■ (32) 

Defining X = —%Y v Uq and supposing that the masses are not degenerate, the solution to 
the previous equation is given by [22,23] 

X = diag(± v /2m 1 M 1 , ±^/2m 2 M 2 , ±^/2m 3 M 3 )/v . (33) 

From this we immediately obtain the constraint (31). Observe that Y v becomes unitary if 
D = 1. However, the present case is not strictly a generalization of the previous one since 
a unitary Y v does not necessarily imply a mass-independent mixing pattern. 

For energy larger than the neutrino mass matrix is fully given by seesaw formula 
(10). The initial condition for m v is given by 

m u(A) = U^rhuUl . (34) 

In the hatted basis J u is proportional to 

Yl% = U D 2 Ul . (35) 

Since Uq is a mass-independent mixing matrix, we should expect that the effect of J v is 
only to change slightly the mass eigenvalues m, but not the mixing angles. In fact, the 
running effect from J v is then proportional to 

2 

(Y^Y U ) T m v (\) + symmetrization = — [/Qdiag(m^Mi, ra\M 2 , rn^M^Ul (36) 

which has exactly the same flavour structure of m v <h\. 

Now we can move to the energy range between and Ms in which the seesaw formula 
is only partial as given in eq. (10). We can exactly proceed in the same way as the previous 

(n) 

case considering first the running effect from J u on k in the hatted basis: 

(«).(«) (n) 

(YjYv) k + symmetrization . 

Explicitly we have 

mj^l k 
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where we have used the unitary condition for Uq. As a result, this contribution is only 
global. 

Now we move to analyze the second term in (10). Observing that in the hatted basis 
we can write 

(n) (n) (n) 

using the unitary condition for C/ , we obtain 

E E( f/ o)^(^o T U(^o) fcm (m) m (t/ t ) mj = J2( U o)kmD 2 m (m) m (U^ mj , 

and the same for the symmetrization part. As we can see, this term also is a global 
contribution. Similarly as in the energy range higher than Ml, also here, the form of m v 
remains invariant and only some of m 8 are slightly shifted. These shifts can be resorbed 
by redefinition of rrii and do not change anyway the mixing angles which are contained in 
Uq and independent from mass eigenvalues. 

Then we arrive to a very general conclusion, in any flavour symmetries with a mass- 
independent mixing pattern, the running effects from J v correct only the neutrino mass 
eigenvalues but not the mixing angles. As in the previous class of models the only flavour- 
dependent RG contribution to ra v is encoded in J e . 

3.1.3 A special case U = iU T BP* and D oc diag(l, 1, —1) 

In this section we consider a special case of Y v = iD Uq in which the expression of Uq is 
enforced by a flavour symmetry based on group. The main feature of this model will be 
presented in the next section together with a more detailed analysis of the running effects. 
Here we need only the constraint on the mixing matrix Uq = iU^bP* and the neutrino 
Yukawa coupling in the hatted basis: 

/ ^/2/3 -1/V6 -1/VE \ 
Y v = yPU^ B 23 = yPi 1/^3 +1/^3 +1/^3 (37) 

V +1/V2 -l/y/2 J 

where y is a positive parameter of order 0(1), P is a diagonal matrix of phases which 
corresponds to the Majorana phases and can be written as 



P = diag(e iQl/2 , e ia2/2 , e ia * /2 ) (31 



and 023 is defined as 



O 



23 




In order to confront (37) with the general expression Y v = iD Uq we observe that 
Y v = yPU^ B 23 U TB U^ B = diagO/,?/, -y)PU$ B . 
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Then we conclude that (37) corresponds to the special case in which D = dia.g(y,y, —y). 
Furthermore, in the model considered in this paper, there is a very simple relation 
between and M; given by mj = v 2 y 2 /2Mi. 

Now we explicitly calculate the RG running from A/ down to A for this special case 
using the approximate analytical expressions given in section 2.1. In the physical basis, 
the light neutrino mass matrix from eq. (3) at the initial energy scale Aj can be recovered 
by imposing the condition in eq. (34): 



m, 



TB 



-U TB Pm u PU± B 








1 








( j 


1 






)+f 1 








l + T ( 






1 





-2 
-2 




(39) 



where m 8 = mie l<H . It is obvious now the meaning of the matrix P as the matrix of the 
Majorana phases of the light neutrinos. It is necessary to specify the kind of neutrino 
mass spectrum: in the Normal Hierarchy (NH) case the light neutrinos are ordered as 
mi < rri2 < ma and the heavy ones as M3 < M2 < Mi; while in the Inverse Hierarchy (IH) 
case they are arranged as 777,3 < m i ^ m 2 and M2 < Mi < M3. 

The general result of the running effects on m u is given by eq. (22) which in our case 
becomes 

m u[X ) = Iu{m T v B + Am' Je) + Am^"') . (40) 

The analytical result for both Ijj and Ami Je ^ (see section 2.1) does not depend on the type 
of the neutrino spectrum, it is sufficient to identify M5, Mm, M^ with the correct hierarchy 
between M!,M 2 ,M 3 . In particular, for the TB mixing pattern, the contribution from J e 
is given by 



777,. 



diag(0, 0, A?) + diag(0, 0, A T )m^ 



TB 








V 



mi 
IT 



777 2 

IT 



777i 






777 2 

IT 



7773 
2 



777,1 


777,2 




\ 


IT 








rrii 


777 2 


7773 




~6~ ~ 


3 


2 




rrii 


2m 2 






3 


3 


- rh 3 


J 



A, 



(41) 



Naturally, the contribution from J v depends on the type of the neutrino spectrum, however 
it can be written in the same form for both the spectra: 



A777, (J " ) 



rn'i 




2777,:, 




where m! i are redefinitions of the light neutrino masses: 



(42) 
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NH case: 



with 



rn'i = mi(p + g) , m' 2 = m 2 (x + q) , mi = m^{x + z) in the SM 

(43) 

fh'i — , wig = 2m 2 x , m 3 = 2m 3 (i + z) in the MSSM 

1 / o 2 ^ 9 2 9 2M M l 

p = -(-3g 2 + AH g 2 + -g 2 )ln — - 

y 16tt 2V y2 KT 1 2 y2; M 2 



1 , „ n . 9 o 9 M 2 



, -3g5 + A H of + -g 2 ) In 

16tt 2V y2 10 Jl 2 y2; M 3 



y 2 1 Mi 

X ~ 32tt 2 M 2 
y2 ln M2 



(44) 



32tt 2 M 



3 



IH case: 



m^ = mi (a: + g) , fh 2 = m 2 (x + 2) , m' 3 = m 3 (p + q) in the SM 

(45) 

m'i = 2mix , mi = 2m 2 (x + z) , ml =0 in the MSSM 



with 



p = - 


1 

16vr 2 


(-3^1 


+ A + 




1 






Q = - 


16vr 2 


(-3^1 


+ A + 


x = - 


y 2 






32vr 2 


Mi 




z = - 


y 2 


, Mi 




32tt 2 


M 2 





9 9 2 .. M 3 

— of + -qi) In 

10 yi 2 y2; Mi 

9 2 9 Mi 
— g? + -g 2 )ln — - 
KT 1 2 y2; M 2 



(46) 



Comparing m^ of eq. (39) with the perturbations Am„ of eqs. (42), we note the 
presence of the same flavour structure for several matrices and in particular, by redefining 
rhi to absorb the terms •m! i it is possible to account for the seesaw contributions from the 
RG running into rn^ B . As a consequence the LO predictions for the TB angles receive 
corrections only from the terms proportional to A T . This result explicitly confirms what 
we outlined in the previous section. 



3.2 RGE effects in the charged lepton sector 

The presence of a term proportional to YjY u in the RG equation for Y e can switch on 
off-diagonal entries in the charged lepton Yukawa matrix Y e . When rotated away, this 
additional contribution introduces a non-trivial U e and consequently corrects the lepton 
mixing matrix Upmns- For a unitary Y v , this correction appears only between the seesaw 
mass scales while, in the general case discussed in section 3.1.2 , it appears already from 
the cutoff Af . 
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In close analogy with the running effects on neutrino mass matrix (40), the full result 
of the running for charged lepton mass matrix can conventionally be written as 



(yX)w = h [(yX)(a } ) + A(y e ty e )] 



(47) 



where I e is an irrelevant global coefficient which can be absorbed by, for example, y T . Now 
we move to the case of TB mixing pattern. In this case, the flavour-dependent corrections 
can be explicitly calculated: 

NH case: 



A(y e t y 



o 








fo 





o\ 




)+c e | 











1 




Vo 





2) 





(48) 



IH case: 



A(y e f y 



^ Vr 





(49) 



where the coefficients are 



16tt 2 3 M 2 ' e 16vr 2 2 M 3 

d = 3CiV * In — a' = -^ V - In ^ 

16tt 2 m SVSY {m z ) ' e 16tt 2 2 M 1 



(50) 



and = -3/2 (1), = 3/2 (3) in the SM (MSSM). Here we observe that the off-diagonal 
contributions to Y^Ye are encoded in a e , b e , a' e and b' e which depend only on the seesaw 
scales Mi . As a result, as we will show in the next section, c e and d e do not affect the 
lepton mixing angles. 



3.3 Full RGE effects on the TB mixing pattern 

In this section, we combine various contributions discussed in previous sections into the 
observable matrix Upmns from which we extract angles and phases at low energy. Since 
we are interested in physical quantities, we eliminate one of the phases of P defined in (38) 
and in particular we express each result in function of = («« — atj)/2, removing a^. The 
corrected mixing angles can be written as 

^K) = e l B + % + ••■ (51) 
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where 9™ = Q, 9\^ = arcsin a/1/3, 9^ 3 = — n/A and fcjj are defined by 



3v/2 V m i - m i 



' 28 



|mi + 777.3!' 

2 2~ 
777.3 _ m l 



1^2 + 777.3]' 
2 2~~ 
777.3 _ m 2 



A r — 3a e — 66 e 



for NH 



777.1 + ^3 



rrin — m 



777,2 + rh 3 \ 



777, — m% 



for IH 



'13 



3V2 



m l — m 3 



777o — m% 



\fhi + m 3 1 

2 2 

mf — 7773 



1^2 + rh 3 \' 

~~~2 ~ 2 
777,o — 777.0 



A T — 3a P 



and the dots stand for sub-leading corrections. In the previous expressions we can clearly 
distinguish the contributions coming from the diagonalization of the corrected TB neutrino 
mass matrix (40) and those from the diagonalization of (47). As it is clear from (50), 
the corrections to the TB mixing from the charged lepton sector is important only for 
hierarchical RH neutrinos and will approach to zero as soon as the spectrum becomes 
degenerate. On the other hand, the corrections from the neutrino sector should be enhanced 
if the light neutrinos are quasi-degenerate and if the tan (3 is large, in the MSSM case. 

The physical Majorana phases are also corrected due to the RG running and we found 
the following results: 

a ij('in x ) — Q ij + daij^r + . . . (52) 

where are the starting values at Af and 

2m 1 m 2 sin(a;i3 - a 23 ) , co s 

0&13 = o 2 2 ( 53 ) 

3 ?77o — 777T 



4 7T7i7T7 2 sin(ai3 - a 23 ) fKA . 
Sa 23 = . (54) 

3 777,2 ~ m l 



At Af, sin^fg 8 is vanishing and as a result the Dirac CP- violating phase is undetermined. 
An alternative is to study the Jarlskog invariants [24] which are well-defined at each energy 
scale: 

JcP — 7^ |^ {{UpMNs)ii(UpMNs)ij(UpMNs)ji(UpMNs)jj} | j (55) 

where i,j G {1,2,3} and i 7^ j. At Af, Jcp is vanishing, while after the RG running it is 
given by 

1 / 777,1 sin ai3 m 2 sin a 23 
Jcp = To 771,3 I — 2 2 2 2 

18 \ 7T7T — 7T7o 777o — 777o 



A r . (56) 
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Two comments are worth. First of all, in the expression for A; 13 , it is easy to recover 
the resulting expression for Jcp as the first term under the square root, a part global 
coefficients. This means that the RG procedure introduce a mixing between the expression 
of the reactor angle and of the Dirac CP-phase. Moreover we can recover the value of the 
Dirac CP-phase directly from eq. (56) and we get the following expression: 



m^mj — 777J) cosai3 — m 2 {m\ — m§) cosa 2 3 — m^{m\ — m%) 

COt OcP = -, o o~V~- 7 2 2~V~- 1" 

mi (rn.2 — m^J sin 0:13 — m2\m\ — m\) sm 023 

3a e (m| — m\){m\ — mf) 
2m3 [r77i(m| — m|) sina 13 — m 2 {m\ — m|) sina^] A T 



(57) 



In the neutrino sector, the RG contributions from the See-Saw terms are present only in 
the resulting mass eigenvalues: 

rrii(x) - + Srrii) + . . . (58) 

where rrii are the starting values at A/ and Srrii, in both the SM and the MSSM and in 
both the NH and IH spectra, are given by 

777i A r „ m'r, 2A T 777o . . , 

8m\ = — ^- , 5m 2 = 2^ - — - , Sm^ = 2 — — A T , 59 

mi 6 7772 3 7773 

with 777^ = 1 771^1, given as in eqs. (43, 45). 



4 The Altarelli-Feruglio (AF) model 

We recall here the main features of the AF model [8], which is based on the flavour group 
Gf = A4 x Z% x U(1)fn'- the spontaneous breaking of A4 is responsible for the TB mixing; 
the cyclic symmetry Z3 prevents the appearance of dangerous couplings and helps keeping 
separated the charged lepton sector and the neutrino one; the U(1)fn [25] provides a nat- 
ural hierarchy among the charged lepton masses. A 4 is the group of the even permutations 
of 4 objects and it has been largely studied in literature [10, 11]: it has 12 elements and 4 
inequivalent irreducible representations, three singlets 1, 1', 1", and a triplet 3. We refer 
to [8] for group details and recall here only the multiplication rules: 



1 ® R = R with R any representation 

1' ® 1' = I" 

I" ® 1" = 1' 

V ®l" =1 

3 ® 3 = 1 © 1' © 1" © 3 S © 3 A • 



(60) 



The TB mixing is achieved through a well-defined symmetry breaking mechanism: A 4 is 
spontaneously broken down to G v = Z 2 in the neutrino sector and to a different subgroup 
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Ge = Z 3 in the charged lepton one. This breaking chain is fundamental in the model, 
because G v and Ge represent the low-energy flavour structures of neutrinos and charged 
leptons, respectively. This mechanism is produced by a set of scalar fields, the flavons, 
which transform only under the flavour group Gf. In table 1, we can see the fermion and 
the scalar content of the model and their transformation properties under Gf. 





i e c fx c t c v c 


H 9 y? T ip s £ 


A, 
Z 3 
U(1) FN 


3 1 1" 1' 3 
UJ u 2 u 2 u 2 u 2 
2 10 


113 3 1 

11 1 UJ 2 UJ 2 

0-10 



Table 1: The transformation properties of the fields under A&, Z% and U(1)fn- 

The specific breaking patter of the symmetry which leads to the TB scheme for leptons 
requires that the flavons develop VEVs with a specific alignment: 

(<Pt) 



(«,0,0) 
c b (u,u,u) 

(61) 

C n U 



A/ 

(<Ps) 
A/ 

<0 

A/ 

(Qfn) _ . 
A/ " 

where c a ^ are complex numbers with absolute value of order one, while u and t are the 
small symmetry breaking parameters of the theory which parametrise the ratio of the 
VEVs of the flavons over the cutoff Af of the theory (they can be taken real through field 
redefinitions). In [8] it has been shown a natural explanation of this misalignment. 

Once defined the transformations of all the fields under Gf, it is possible to write down 
the Yukawa interactions: at the LO they read 

Cf° = ^9 2 e c H^ip T e) + ^efx c H^^ T £y + ^T c HHip T £)'' + h.c. (62) 

C L v ° = y(v cT H^,) + x a Z{v cT v c )+x b { Vs v cT v c ) + h.c. . (63) 

The notation (...), (•••)' and (•••)" refers to the contractions in 1, V and 1", respectively. 
When the flavons develop VEVs in agreement with eq. (61) and after the electroweak 
symmetry breaking, the LO mass matrix of charged leptons takes the following form: 

t 2 
m t = ^\ y„t 
y T 
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with y e , Un and y T being complex numbers with absolute values of order one. As one can 
see the relative hierarchy among the charged lepton masses is given by the parameter t: 
when 

t « 0.05 (65) 

then the mass hierarchy is in agreement with the experimental measurements. As we will 
see in the following sections, the model admits a well defined range for the parameter u 
which can approximatively be set to 

0-003 <u< 0.05 . (66) 

In the neutrino sector, the Dirac and the Majorana mass matrices, at the LO, are given 

by 

/ 1 \ / a + 26/3 -6/3 -6/3 \ 

m D = ^ I 1 M R oc -6/3 26/3 a - 6/3 , (67) 

^ 2 \ 1 / \ -6/3 a - 6/3 26/3 / 

where a = 2x a c a u and 6 = 2xbCbU. The complex symmetric matrix Mr is diagonalised by 
the transformation 

M R = U t r MrUr , (68) 
where Mr is a diagonal matrix with real and positive entries, given by 

Mr = diag(M 1; M 2 , M 3 ) = diag(|a + b\,\a\, \ - a + b\) , (69) 

while the unitary matrix Ur can be written as, 

U R = U TB P , (70) 

where P is the diagonal matrix of the Majorana phases already defined in eq. (38), with 
«i = — arg(a + 6), a 2 = —arg(a), a 3 = — arg(— a + b). After the electroweak symmetry 
breaking, the mass matrix for the light neutrinos is recovered from the well known type I 
seesaw formula 2 

m v = -m T D M^m D = - V -Jj-M~ l (71) 

where the last passage is possible considering that M^tud = moM^ 1 . From eq. (68), 
U R M R l U R = diag(Mf \ M 2 _1 , M^ 1 ) and clS db result the light neutrino mass matrix can be 
diagonalised by 

m„ = U T v m v U v , (72) 

where U u = iU R = iUtbP* and the diagonal matrix m„ has real and positive entries written 

as 5 

v 2 y 2 

mi = il ' (73) 



5 Note that we can absorb the phase of y in P without loss of generality: it is a global non-observable 
phase. 
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When the SUSY context is considered Gf accounts for an additional term, a continuous 
R symmetry U(1)r which simplifies the constructions of the scalar potential: under this 
symmetry, the matter superfields are single-valued, while the scalar ones are neutral. It is 
easy to extend eqs. (62,63) in the supersymmetric case: two Higgs doublets hd, u , invariant 
under A±, substitute H and H, respectively; the Lagrangian Cf° is identified to the LO 
charge lepton superpotential w\° and is identified to the LO neutrino superpotential 
w^°. While t is still equal to 0.05 in order to have a correct charged lepton mass hierarchy, 
the range for u slightly changes: 

0-007 < u < 0.05 . (74) 



4.1 The light neutrino mass spectrum 

Following [14,19], we summarize the results for the light neutrino mass spectrum, which 
will be useful to discuss the effects of the running on the low-energy observables. We note 
that the following analysis is valid in the SM as well as in the MSSM. The light neutrino 
masses are directly linked to the heavy neutrino masses through eq. (73) and they can be 
expressed in terms of only three independent parameters. It is possible to choose these 
parameters in order to simply the analysis: \a\ = \M 2 \ = v 2 \y\ 2 /(2\m 2 \) (in the MSSM we 
replace v with v u ), p and A, where p and A are defined as 

b 



a 



pe lA , (75) 



where A is defined in the range [0, 2tt]. From the experimental side only the squared mass 
differences have been measured: for the NH (IH) they are [5] 

Am s 2 ol = m\ — m\ 

(76) 

A ™atm = \m\-m\[m\)\ . 

As a result the spectrum is not fully determined and indeed A is still a free parameter. 
We can bound this parameter, requiring | cos A| < 1. In order to get analytical relations 
for p and cos A, we calculate the following mass ratios: 

2 

= l±2pcosA + p 2 . (77) 



m i(3) 



It is then easy to express p and cos A as a function of the neutrino masses: 




2 2 
77*2 

2 2 

cos A = 1 1 . (7£ 




20 



It is interesting to note that this expression holds for both the types of spectra. Using 
eqs. (76), it is possible to express cos A in function of only the lightest neutrino mass 
and, imposing the constraint | cos A| < 1, the following ranges can be derived: taking the 
central values of Am s 2 ol and Am^ tm 

NH : 4.46 meV < m 1 < 5.91 meV 

(79) 

IH : 17.1 meV < m 3 . 

For the NH, mi spans in a narrow range of values, which corresponds to values of A close 
to zero. On the other hand, for the IH, m 3 is bounded only from below and the minimum 
is achieved when A is close to ±7r. Furthermore, cos A is restricted only to negative values 
(see [14,19] for further details). 




io- 3 i(r 2 ltr 1 i 10" 3 10" 2 lo- 



rn i (eV) m-i (eV) 

Figure 1: Plots of the ratios of the heavy RH neutrino masses as a function of the lightest LH neutrino 
masses. On the left the NH case and on the right the IH one. The green opaque areas refer to the allowed 
range for mu$\ as in eq. (79). The vertical black lines correspond to the future sensitivity of 0.2 eV of 
KATRIN experiment. 

From eq. (73) it is possible to describe the LO spectrum of the RH neutrinos in function 
of a unique single parameter, which is the lightest LH neutrino mass. In all the allowed 
range for m 1 ( 3 ), the order of magnitude of the RH neutrino masses is 10 14 ^" 15 GeV. In fig. 
(1) we show explicitly the ratios of the RH neutrino masses for NH and the IH cases, on 
the left and on the right respectively. The ratios are well defined for the NH, thanks to 
the narrow allowed range for m\\ M1/M3 ~ 11 and M 2 /M 3 ~ 5. In the case of the IH, 
the ratio M1/M2 is fixed at 1 while M3/M2 varies from about 3 to 1, going from the lower 
bound of m 3 up to the KATRIN sensitivity. 

These results are valid only at LO and some deviations are expected with the introduc- 
tion of the higher order terms, that is illustrated in the following section. The result of a 
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direct computation shows that for the NH spectrum the corrections leave approximatively 
unaffected eqs. (79); this is true for the IH case too, apart when the neutrino masses reach 
values at about 0.1 eV for which the deviations become significant. 



We only commented on the neutrino masses, both light and heavy, and we saw that it 
is possible to express all of them in function of the lightest LH neutrino mass. The same is 
true for the phases too: since in the TB mixing the reactor angle is vanishing, the Dirac CP 
phase is undetermined at LO; on the contrary the Majorana phases are well defined and 
they can be expressed through p and A. Since we are interested in physical observables, 
we report only phases differences, = (ctj — <x,)/2: 



• ro \ psinA 
sm(2a 23 ) = 



sin (2a 



13 J 



y/l - 2pcosA + p 2 

2psin A 
vV - l) 2 + 4p 2 sin 2 A 



(80) 



It will be useful for the following discussion to show also sin(2«i2), which enters in the RG 
evolution of the physical Majorana phases: 

• /o \ psinA 
sm(2a 12 ) = - 



1 + 2pcos A + p 2 



4.2 The next-to-leading order (NLO) contributions 

Another important implication of the spontaneously broken flavour symmetry is that the 
leading order predictions are always subjected to corrections due to higher-dimensional 
operators. The latter are suppressed by additional powers of the cutoff Aj and can be 
organized in a suitable double power expansion in u and t. 

At the NLO there are may additional terms which can be added in the Lagrangian. 
Since <pt is the only scalar field which is neutral under the abelian part of the flavour 
symmetry, all the NLO terms contain the terms already present in the LO Lagrangian 
multiplied by an additional ipx/Af . In addition to these terms, there are also corrections 
to the leading vacuum alignment in eq. (61): 



A/ 
A/ 

<£> 

A/ 



(w,0,0) + {c x u 2 , c 2 u 2 , c 3 u 2 ) 

c b (u, u, u) + (c 4 m 2 , c 5 m 2 , c e u 2 ) (82) 



c a u + c 7 u 2 



where Cj are complex numbers with absolute value of order one. Note that in the MSSM, 
the model predicts c 2 = C3. Here we will not perform a detailed analysis for NLO operators 



22 



and the origin of (82), see [8,19] for a relative study. As a result, all the quantities relevant 
in this paper receive a NLO correction at the 0(u) level: 



Y e + 5Y e , Y u + 8Y V , M + 5M . 

These corrections should modify the resulting m v and the mixing angles in addition to the 
RG contributions. It is not difficult to see that deviations from TB induced only by NLO 
corrections are also of the same level [8, 19]: 

sin 2 9 23 = - + 0{u), sin 2 On = \ + 0(u), sin 9 13 = 0{u). (83) 

Since the solar mixing angle is, at present, the most precisely known, we require that its 
value remains at 3er range [5]. This requirement results in an upper bound on u of about 
0.05. On the other hand we have from eq. (64) the following relation: 



u = ~ 0.01 A in the SM 

IVrl v \y T \ 

tan8\/2m T _ tan/3 . , 

u ~ — j — j-— « 0.01— ; — V- in the MSSM 

2/t v \y T \ 



where for the r lepton we have used its pole mass m T = (1776.84 ± 0.17) MeV [18]. 
Requesting \y T \ < 3 we find a lower limit for u of about 0.003 in the SM case; in the 
MSSM, the same requirements provides a lower bound close to the upper bound 0.05 for 
tan f3 = 15, whereas for tan (3 = 2 it is u > 0.007. From now on, we will choose the maximal 
range of u as 

0.003 < u < 0.05 (85) 
for the SM context, while for the MSSM one we take 

0.007 < u < 0.05 , (86) 

which shrinks when tan/3 is increased from 2 to 15. 

The NLO terms affect also the previous results for the Dirac and the Majorana phases. 
All the new parameters which perturb the LO results are complex and therefore they 
introduce corrections to the phases of the PMNS matrix. Due to the large amount of such 
a parameters, we expect large deviations from the LO values. 



5 RG effects in the Altarelli-Feruglio model 

In this section we will apply the analysis of RG running effects on the lepton mixing angles 
to the AF model. In order to perform such a study, it is important to verify the initial 
assumptions made in section 3.3, in particular, we see that eq. (37) exactly corresponds to 
the one implied by the AF model, when moving to the physical basis. On the other side, the 
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presence of flavon fields has a relevant impact on the results of the analysis. In the unbroken 
phase, flavons are active fields and should modify the RG equations. Although the study 
of the relevant Feynman diagrams goes beyond the aim of this work, what follows can be 
easily proved. Since the only source of the A± breaking is the VEVs of the flavons, any 
flavour structure is preserved above the corresponding energy scale, whatever interactions 
are present. In particular, the lagrangian (62, 63) contains all possible LO terms, given 
the group assignments, and its invariance under A 4 is maintained moving downward to the 
scale (ip), where significant changes in the flavour structure can appear. From eqs. (67) and 
(69), we deduce that (<p) ~ Mi and as a result in the AF model A T must be proportional 
to ]n(((p)/X) and not to ln(Aj/A). 

We will separately discuss the evolution of angles and phases for both type of hierarchy. 
In the following, the results will be shown for the SM and for the MSSM with tan j3 = 15 
apart where explicitly indicated otherwise. Without loss of generality, we choose y = 1 for 
our numerical analysis. We also set (if) = 10 15 . The spectrum spans the range obtained 
in (79). 

5.1 Running of the angles 

Since we are interested in deviations of the corrected mixing angles from the TB predictions 
and in confronting them with experimental values, it is convenient to relate the coefficients 
kij defined in section (3.3) with physical observables. Keeping in mind that \kij\ <C 1 and 
that we start from a TB mixing matrix, it follows that 

1 2 v^2 

sin 6>i3 ~ k 13 , cos 26 23 ~ 2k 23 , sin 2 6 12 - - ~ -3-^12 • ( 87 ) 

The corrections to the TB mixing angles as functions of mi(m 3 ) in the NH (IH) case are 
shown in figure 2. 

We begin with the case of NH. Since the dependence of the corrected mixing angles 
from A r is the same, SM corrections are generally expected to be smaller than those in 
MSSM. However, from figure 2 we see that, in NH, there is not a large split between the two 
curves for SM and MSSM respectively. This fact suggests a dominant contribution coming 
from the charged lepton sector. For the atmospheric and reactor angles, the deviation from 
the TB prediction lays roughly one order of magnitude below the la limit. In particular, 
RG effects on sin 8 r3 are even smaller than the NLO contributions analyzed in section 4.2 
which are of 0(u), without cancellations. On the other hand, since the experimental value 
of the solar angle is better measured than the other two, the running effects become more 
important in this case. Indeed, RG correction to the TB solar angle evades the la limit as 
it can be clearly seen in figure 2. Anyway, we observe that for both the atmospheric and 
solar angles, the running contribution is of the same order as the contribution from NLO 
operators. 

Now we move to analyze the case of IH. In this case, since the neutrino spectrum 
predicted by the AF model is almost degenerate, the contribution from the charged lepton 
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Figure 2: Corrections to the TB mixing angles as functions of 7711(7713) for the NH (IH) are shown. On the 
left column the plots refer to the NH spectrum, while on the right column the IH case is reported. The 
plots show the MSSM case with tan/3 = 15 (solid blue) and the SM case (black dashed), compared to the 
current la and 3cr limits (dashed red). For the specific case of | sin 2 612 — 1/3| for the IH, two values for 
tan/3 — 2, 5, 9 are considered (solid blue, purple and green). All the plots span in a range for 7711(7713) 
which is given by eq. (79) or by the KATRIN bound. 



sector in eqs. (49) is subdominant. As a consequence the information which distinguishes 
the SM case from the MSSM one is mainly dictated by A r defined in eq. (25). As a result 
the running effects in the MSSM are always larger than in the SM and for large tan (3 they 
are potentially dangerous. The curves corresponding to the atmospheric and reactor angles 
do not go above the 3a and la windows respectively. However, the deviation from 
presents a more interesting situation. For example, for tan/3 > 10, the RG effects push 
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the value of the solar angle beyond the 3a limit for entire spectrum. For lower values of 
tan /3, the model is within the 3a limit only for a part of the spectrum where the neutrinos 
are less degenerate. Confronting with the running effects, in the IH case, the contribution 
from NLO operators in the AF model is under control. 

5.2 Running of the phases 

Majorana phases are affected by RG running effects too. Since there is not experimental 
information on Majorana phases available in this moment we will simply show their values 
at low energy, eventually comparing them with the prediction in the AF model. We stress 
again that they are completely determined by only one parameter, the mass of the lightest 
neutrino, mi for NH and for IH. 




0.0046 0.0050 0.0054 0.0058 0.02 0.06 0.10 0.14 0.18 

mi (eV) m 3 (eV) 

Figure 3: Majorana phases ai3 and «23 as functions of the lightest LH neutrino masses. For the NH (left 
panel) the corresponding curves at low and high energies are undistinguishable. For the IH (right panel) 
the curves refers to low energy values in MSSM with tan/3 =15 (solid blue or red) and the AF prediction 
at A/ (dashed blue or red). 

In the case of NH, Majorana phases are essentially not corrected by RG effects. This 
feature is due to the fact that 5a±3 and 5a 2 3 of eqs. (53, 54) are proportional to sin(ai 3 — a 23 ) 
which is close to zero, as we can see looking at the left panel of figure 3. In the case of IH, 
MSSM RG effects always increase the values of phases when moving from high energy to 
low energy and they are maximized for tan/3 = 15, especially when the neutrino spectrum 
becomes degenerate. On the contrary, in the SM context, the low energy curves cannot be 
distinguished from the high energy ones. 

As described in section (3.3), a definite Dirac CP violating phase 5c p arises from 
running effects even if, in the presence of a TB mixing pattern, it is undetermined in the 
beginning. Although the final Dirac phase can be large, Jarlskog invariant, which measures 
an observable CP violation, remains small because of the smallness of 6* 13 . We remember 
that these results are valid both for the SM and for MSSM. 



26 



100 

80 
60 
40 
20 




^ 10 



0.0046 




0.0046 0.0050 0.0054 
m\ (eV) 



0.0058 



0.02 



0.06 



0.10 0.14 

in?, (eV) 



0.0050 0.0054 
mi (eV) 



0.0058 



0.02 0.06 



0.10 0.14 
;«3 (eV) 



0.18 




0.18 



Figure 4: Dirac CP phase Sep and Jaiiskog invariant Jcp as functions of the lightest LH neutrino masses 
for the NH (left column) and the IH (right column), in the MSSM with tan/3 = 15. 



6 Conclusions and outlook 

The flavour sector is poorly known within the SM or the MSSM and the masses and mix- 
ings are considered as free parameters that can be adjusted to agree with experiments. 
This lack can be improved by adding to the SM or the MSSM gauge group appropriate 
flavour symmetries such that Yukawa couplings can be understood in a more fundamental 
way. Even though the mechanism that generates the fermion masses is not yet completely 
understood, a great effort has been made in the last years, especially after the discovery of 
neutrino masses which are significantly smaller than those of charged fermions. Further- 
more, the leptonic mixing pattern is also very different from Vckm because it contains a 
nearly maximal atmospheric angle #23 and a solar angle very approximate to the TB pre- 
diction. The flavour structure of the neutrino mass matrix can be nicely explained by (type 
I) seesaw mechanism implemented with flavour symmetries. In order to naturally describe 
lepton masses and mixings, however, the flavour group must be (spontaneously) broken. 
Then we can expect that the LO prediction of a flavour symmetry is always subject by 
subleading corrections characterized by small symmetry breaking parameters (t and u in 
our case). Then, in a consistent flavour model building, these subleading corrections must 
be safely under control. 
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In any flavour model, there is also an other independent correction to the LO predictions 
which is due to the RG evolution of parameters. In this paper we have studied these running 
effects on neutrino mixing patterns. In seesaw models, the running contribution from the 
neutrino Yukawa coupling Y v , encoded in J„ , is generally dominant at energy above the 
seesaw threshold. However, this effect, which in general introduces appreciable deviations 
from the LO mixing patterns, does not affect the mixing angles, under specific conditions: 
in the first part of the paper, we have analyzed two classes of models in which this indeed 
happens. The first class is characterized by a unitary Y v . It is the case, for example, 
when the RH singlet neutrinos or the charged leptons being in a irreducible representation 
of the flavour group. The second class is the so called mass-independent mixing pattern, 
in which, in particular, the effect of J„ can be absorbed by a small shift on neutrino 
mass eigenvalues leaving mixing angles unchanged. The widely studied TB mixing pattern 
belongs, for example, to this second class of models. 

In the second part of the paper, we focused on a special realization of the general class 
of flavour models studied in the first part. We were interested in the AF model for TB 
mixing where the flavour symmetry group is given by A4 x Z% x {7(1)fn- The aim is to 
analyze the RG effects on the TB mixing pattern in addition to the NLO corrections already 
present in this model and to confront them with experimental values. The analysis has 
been performed both in the SM and MSSM and for both neutrino spectra. We found that 
for NH light neutrinos, the dominant running contribution comes from the charged lepton 
sector which weakly depends on both tan/3 and mass degeneracy. As a result, for this 
type of spectrum, the tribimaximal prediction is stable under RG evolution. Moreover, the 
running contribution is of the same order or smaller with respect to the contribution from 
NLO operators. On the other hand, in the case of IH, the deviation of the solar angle from 
its TB value can be larger than the NLO contribution and, in particular, for tan f3 > 10 an 
IH spectrum is strongly disfavored. In the end, observe that for both spectra, the reactor 
angle #13 does not receive appreciable deviations from zero (at a level < u). 

The effects of RG running can be manifested also in other phenomena which are not 
directly related to the neutrino properties, such as Lepton Flavour Violating (LFV) tran- 
sitions, leptogenesis, etc. For example, in [19,22,23,26], it has been pointed out that, in 
the limit of an exact A4 symmetry, all CP violating asymmetries vanish. Then when the 
flavour symmetry is spontaneously broken, NLO corrections become important in gener- 
ating the desired leptogenesis. On the other hand, also the RG evolution can introduce 
symmetry breaking effects which can be in principle dominant over NLO contributions. As 
a result, the estimate of generated leptogenesis can be quite different from those obtained 
in [19,22,23,26] since they do not take into account the RG effects. Another important 
consequence of RG running is the generation of off diagonal terms in the soft SUSY break- 
ing mass matrices, contributing to LFV rare processes. In a series of papers [27-29] it 
has been studied the impact of using flavour symmetries in order to explain the measured 
bounds on rare decays. In [27] it has already been shown that, below the seesaw scales, 
the running effect is negligible with respect to that originated in the corresponding flavour 
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theory. However between the seesaw scales the threshold effects are important and the two 
contributions to LFVs can be comparable. All these issues are very interesting and are 
subject for a further investigation. 
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Appendix A : Renormalisation group equations 



In order to calculate the RG evolution of the light neutrino mass matrix, the RGEs for 
all the parameters of the theory have to be solved simultaneously. We use the notation 
defined in the text, where a superscript (n) denotes a quantity between the nth and the 
(n + l)th mass threshold. When all the RH neutrinos are integrated out, the RGEs can 
be recovered by setting the neutrino Yukawa coupling Y v to zero, while in the full theory 

above the highest seesaw scale, the superscript (n) has to be omitted. 

(«•) (™) („) 

In the SM extended by singlet neutrinos, the RGEs for Y e , Y u , M, k, Yd, and Y u are 
given by 



(») f 3 , 3 (")>) 

7T 2 F R = Y R \ —Y^Y e — —YjY v + Tr 



(n)» 



3 YjY u + 3 YjY d + YjY u + YjYe 



4 0i A 9l 



An) (») r 3 («)(») 3 

\^y v = y t J -yX - -y e ty e 



Tr 



(«)>) 



3 YlY u + 3 YjY d + Y}Y V + Y}Y t 
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(n) 



16vr 2 K 



(n) In) 
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P (n) l(n) 

ft H K 

2 



(«) (n) 



+ 



+2Tr 



(«)» 



3 y^y u + 3 y;y d + yJy v + y e ty 



- 3g$K + X h k , 



3 yJYu + 3 yjy + y)y v + y;y e 



{3 3 j. r 4. («) (») 

-y n ty u - -yjy, + Tr [3 y u ty u + 3 yjy, + yJy v + yjy 



1 2 
4^ 
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17 9 
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(n) 



1Qtt z X h 



QX 2 h - 3A„ ( S^ 2 + ~<? 2 ) + 3g\ + ? + s 2 



(n).(rt) 



3 yjy u + 3 y]y d + yjy, + y e ty 



+4A# Tr 
-8Tr 

We use the convention that the Higgs self-interaction term in the Lagrangian is —\h{H^H) 2 /4:. 



. , (n) (n) (n) (n) 

3 yjy tt yjy u + 3 y\y a yJy, + yJy u yX + rfr y e ty e 
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In the MSSM context the 1-loop RGEs for the same quantities are given by 



d i n ) 



167r 2 -F e = Y e { 3Y}Y e + YX + Tr 3Y}Y d + Y}Y e - § g\ - 3g 
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